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Abstract 

We present new measurements of CP-violation parameters in B° — > (j)K , K + K~ Kg, fo(980)K s , 
rj'K s , ujK s , and K s ir° decays based on a sample of 275 x 10 6 BB pairs collected at the T(45) 
resonance with the Belle detector at the KEKB energy-asymmetric e + e~ collider. One neutral B 
meson is fully reconstructed in one of the specified decay channels, and the flavor of the accompa- 
nying B meson is identified from its decay products. CP- violation parameters for each of the decay 
modes are obtained from the asymmetries in the distributions of the proper-time intervals between 
the two B decays. The combined result for the B° -► (j)K°, K+K~K S , f (980)K^, r)'K Q s , loK q s , 
Kgir and previously reported KgKgKg decays differs from the SM expectation by 2.4 standard 
deviations. 

PACS numbers: 11.30.Er, 12.15.Hh, 13.25.Hw 



3 



I. INTRODUCTION 



The flavor- changing b — > s transition proceeds through loop penguin diagrams. Such 
loop diagrams play an important role in testing the standard model (SM) and new physics 
because particles beyond the SM can contribute via additional loop diagrams. CP violation 
in the b — > s transition is especially sensitive to physics at a very high-energy scale [1]. 
Theoretical studies indicate that large deviations from the SM expectations are allowed for 
time-dependent CP asymmetries in B° meson decays [2]. Experimental investigations have 
recently been launched at the two B factories, each of which has produced more than 10 8 
BB pairs. Our previous measurement of the CP- violating asymmetry in B° — > <f)Kg [3], 
which is dominated by the b — > sss transition, yielded a value that differs from the SM 
expectation by 3.5 standard deviations [4]. Measurements with a larger data sample are 
required to elucidate this difference. It is also essential to examine additional modes that 
may be sensitive to the same b — > s penguin amplitude. In this spirit, experimental results 
using the decay modes B° -> <f)K° L , K + K~K° S , f (980)K° s , rfK° s , and K° s tc° have already 
been reported [4, 5]. 

In the SM, CP violation arises from a single irreducible phase, the Kobayashi-Maskawa 
(KM) phase [6], in the weak-interaction quark- mixing matrix. In particular, the SM predicts 
CP asymmetries in the time- dependent rates for B° and B° decays to a common CP eigen- 
state fcp [7]. In the decay chain T(AS) — > B°B° — > fcpftag, where one of the B mesons 
decays at time t C p to a final state fcp and the other decays at time t tag to a final state / tag 
that distinguishes between B° and B°, the decay rate has a time dependence given by 

e -\At\/r B o f n 

V(At) = \l + q- Ssm(Am d At) + Acos(Am d At) \. (1) 

Here S and A are CP- violation parameters, t b ° is the B° lifetime, Am d is the mass difference 
between the two B° mass eigenstates, At = tcp — t tag , and the 6-flavor charge q — +1 
(— 1) when the tagging B meson is a B° (B°). To a good approximation, the SM predicts 
S = — £/ sin 20i, where £/ = +1(— 1) corresponds to CP-even (-odd) final states, and ^4 = 
for both b — > ccs and b — > sqq transitions. Recent measurements of time- dependent CP 
asymmetries in B° — > J/ipK^ and related decay modes by Belle [8, 9] and BaBar [10], which 
are governed by the b — > ccs transition, have already determined sin 20i rather precisely; the 
present world average value is sin20x = +0.726 ± 0.037 [11]. This serves as a firm reference 
point for the SM. 

Belle's previous measurements for B° — > (f)Kg, K + K~Kg and rf'Kg were based on a 140 
fb^ 1 data sample (DS-I) containing 152 x 10 6 BB pairs. In this report, we describe improved 
measurements incorporating an additional 113 fb -1 data sample that contains 123 x 10 6 BB 
pairs (DS-II) for a total of 275 x 10 6 BB pairs. 

While (j)Kg and rj'Kg final states are CP eigenstates with £j = —1, the K + K~K° S final 
state is in general a mixture of both £j = +1 and —1. Excluding K + K~ pairs that are 
consistent with a <fi — > K + K~ decay from the B° — > K + K~Kg sample, we find that the 
K + K~Kg state is primarily £j = +1; a measurement of the = +1 fraction using the 
isospin relation [12] with a 253 fb" 1 data sample gives /+ = 0.83 ± O.lO(stat) ± 0.04(syst), 
which is consistent with the previous result [4]. The SM expectation for this mode is S = 
-(2/+ -1) sin 20i. 

We include additional <pKg and rj'Kg subdecay modes that were not used in the previous 
analysis. We also describe new measurements of CP asymmetries for the following CP- 
eigenstate B° decay modes: B° -> (\>K\ and f (980)K° s for ^ = +1; B° -> uK° s and K° s tt 
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for £f = — 1. The decays B° — > <fiK s and (pK^ are combined in this analysis by redefining 
S as —£fS to take the opposite CP eigenvalues into account, and are collectively called 
U B° — > <f)K°" . The CP asymmetries for the decay B° — > w-ffg are measured for the first 
time. 

At the KEKB energy- asymmetric e + e~ (3.5 on 8.0 GeV) collider [13], the T(4S) is pro- 
duced with a Lorentz boost of /?7 = 0.425 nearly along the electron beamline (z). Since 
the B° and U° mesons are approximately at rest in the T(4S) center-of-mass system (cms), 
At can be determined from the displacement in z between the fcp and / tag decay vertices: 
At ~ (z CP - z tag )/(P 1C ) = Az/i^c). 

The Belle detector is a large-solid-angle magnetic spectrometer that consists of a silicon 
vertex detector (SVD), a 50-layer central drift chamber (CDC), an array of aerogel threshold 
Cherenkov counters (ACC), a barrel-like arrangement of time-of-flight scintillation counters 
(TOF), and an electromagnetic calorimeter comprised of CsI(Tl) crystals (ECL) located 
inside a superconducting solenoid coil that provides a 1.5 T magnetic field. An iron flux- 
return located outside of the coil is instrumented to detect K\ mesons and to identify muons 
(KLM). The detector is described in detail elsewhere [14]. Two inner detector configurations 
were used. A 2.0 cm radius beampipe and a 3-layer silicon vertex detector (SVD-I) were 
used for DS-I, while a 1.5 cm radius beampipe, a 4-layer silicon detector (SVD-II) and a 
small-cell inner drift chamber were used for DS-II [15]. 



II. EVENT SELECTION, FLAVOR TAGGING AND VERTEX RECONSTRUC- 
TION 

A. Overview 

We reconstruct the following B° decay modes to measure CP asymmetries: B° — > (f)K s , 
<f)K° L , K + K~K S , f (980)K° s , rfK%, uK° s , and K s tt . We exclude K + K~ pairs that are 
consistent with a <fi — > K + K~ decay from the B° — > K + K~K S sample. The intermediate 
meson states are reconstructed from the following decays: it — > 77, K s — > it + it~~ (also 7r°7r° 
for the (pKg decay), i] — > 77 or 7r + 7r~7r°, p° — > 7r + 7r~, uj — > n + n^7Y°, rj' — > p°7 or r]n + n^, 
/o(980) -> 7T+7T-, and <p -> 

Among the decay chains listed above, 5° — > 0^ (ifg — > B° — > K + K~Kg, 

B° — > V-^s I 7 ?' ~" *■ P°7)> an d ~" *■ ^'^5 I 1 ?' ~ ^7r + 7r _ , 77 — > 77) decays were used in the 
previous analysis [4]. The selection criteria for these decays remain the same. For the newly 
included B° decay modes, identification of photons, neutral and charged kaons, and neutral 
and charged pions is based on the procedure used previously. However, the selection criteria 
for each B° decay mode were optimized individually and are thus different from one another. 



B. B° -» (j)K% and K+K~K% 

Charged tracks reconstructed with the CDC for kaon and pion candidates except for 
those from Kg — > tt + tt~ decays are required to originate from the interaction point (IP). We 
distinguish charged kaons from pions based on a kaon (pion) likelihood Ck(tt) derived from 
the TOF, ACC and dE/dx measurements in the CDC. 

Pairs of oppositely charged tracks that have an invariant mass within 0.030 GeV/c 2 of 
the nominal K° s mass are used to reconstruct Kg — > it + it~ decays. The distance of closest 
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approach of the candidate charged tracks to the IP in the plane perpendicular to z axis is 
required to be larger than 0.02 cm for high momentum (> 1.5 GeV/c) K$ candidates and 
0.03 cm for those with momentum less than 1.5 GeV/c. The 7r + 7r~ vertex is required to 
be displaced from the IP by a minimum transverse distance of 0.22 cm for high-momentum 
candidates and 0.08 cm for the remaining candidates. The mismatch in the z direction at the 
Kg vertex point for the 7r + 7r~ tracks must be less than 2.4 cm for high-momentum candidates 
and 1.8 cm for the remaining candidates. The direction of the pion pair momentum must 
also agree with the direction of the vertex point from the IP to within 0.03 rad for high- 
momentum candidates, and to within 0.1 rad for the remaining candidates. The resolution 
of the reconstructed K s mass is 0.003 GeV/c 2 . 

Photons are identified as isolated ECL clusters that are not matched to any charged 
track. To select Kg — > 7r°7r° decays, we reconstruct ir° candidates from pairs of photons 
with > 0.05 GeV, where E 1 is the photon energy measured with the ECL. Photon pairs 
with an invariant mass between 0.08 and 0.15 GeV/c 2 and a momentum above 0.1 GeV/c are 
used as ir° candidates. Initially, the 7r° decay vertex is assumed to be the IP. The asymmetric 
mass window is used to take into account the lower tail of the mass distribution due to the 
distance between the IP and the true ir° vertex. Candidate K s — > 7r°7r° decays are required 
to have an invariant mass between 0.47 GeV/c 2 and 0.52 GeV/c 2 , where we perform a fit 
with constraints on the K s vertex and the 7r° masses to improve the 7r°7r° invariant mass 
resolution. We also require that the distance between the IP and the reconstructed K s decay 
vertex be larger than —10 cm, where the positive direction is defined by the K s momentum. 

Candidate — > K + K~ decays are required to have an invariant mass that is within 
0.01 GeV/c 2 of the nominal meson mass. Since the meson selection is effective 
in reducing background events, we impose only minimal kaon-identification requirements; 
Tlx/n = £>k/(£>k + C-k) > 0.1 is required, where the kaon likelihood ratio TZr/tt has values 
between (likely to be a pion) and 1 (likely to be a kaon). We use a more stringent kaon- 
identification requirement, TZx/n > 0.6, to select non-resonant K + K~ candidates for the 
decay B° — > K + K~ K s . The K + K~ candidates for B° — > K + K~K S are selected by reject- 
ing K + K~ pairs with an invariant mass within 0.015 GeV/c 2 of the nominal meson mass, 
reducing the contribution to a negligible level. To remove Xco ~ * K + K~, J/ip — > K + K~ 
and D° — > K + K~ decays, K + K~ pairs with an invariant mass within 0.015 GeV/c 2 of the 
nominal masses of Xco and J/ip or within 0.01 GeV/c 2 of the nominal D° mass are rejected. 
D + — > K° S K + decays are also removed by rejecting K S K + pairs with an invariant mass 
within 0.01 GeV/c 2 of the nominal D + mass. 

For reconstructed B — > fcp candidates, we identify B meson decays using the en- 
ergy difference AE = E c ™ — E™* m and the beam-energy constrained mass M bc = 

\J(Eb™m) 2 ~ (PB ns ) 2 > where E^ m is the beam energy in the cms, and E c ™ and p c ^ s are 
the cms energy and momentum of the reconstructed B candidate, respectively. The reso- 
lution of Mbc is about 0.003 GeV/c 2 . Because of the smallness of p ™, the Mb c resolution 
is dominated by the beam-energy spread which is common to all decay modes. The resolu- 
tion of AE depends on the reconstructed decay mode. The AE resolution is 0.013 GeV 
for (j)K s {K° s -> 7T+7T-) and K+K~K S . The AE distribution for <pK s {K° s -> 7r°7r°) 
has a tail toward lower AE due to 7 energy leaking in the ECL. The AE resolution for 
4>K S (K s — > 7r°7r°) is 0.058 GeV for the main component and the width of the tail com- 
ponent is about 0.14 GeV. The B meson signal region is defined as \AE\ < 0.06 GeV for 
B° -> (j)K s (K s -> 7T+7T-), -0.15 GeV < AE < 0.1 GeV for B° -> 07^ (K° s -> n°n ), 
\AE\ < 0.04 GeV for B° -> K + K~K° S , and 5.27 GeV/c 2 < M bc < 5.29 GeV/c 2 for all 
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decays. 

The dominant background to the B° — > <pK s decay comes from e+e" — > ml, rfrf, ss, 
or cc continuum events. Since these tend to be jet-like, while the signal events tend to 
be spherical, we use a set of variables that characterize the event topology to distinguish 
between the two. We combine S±, 9t and modified Fox- Wolfram moments [16] into a Fisher 
discriminant J 7 , where S± is the scalar sum of the transverse momenta of particles other than 
the reconstructed B candidate outside a 45° cone around the candidate <p meson direction 
(the thrust axis of the B candidate for K + K~K S decays) divided by the scalar sum of their 
total momenta, and Ot is the angle between the thrust axis of the B candidate and that 
of the other particles in the cms. We also use the angle of the reconstructed B candidate 
with respect to the beam direction in the cms (9b), and the helicity angle 9h defined as 
the angle between the B meson momentum and the daughter K + momentum in the (p 
meson rest frame. We combine J 7 , cos9 B and cos9 H into a signal [background] likelihood 
variable, which is defined as £ S i g [bk g ] = £ S i g [bk g ] (F) x £si g [bk g ](cos#e) x £ s i g [bk g ] (cos 9 H ) ■ For 
K + K~K S decays, T and cos9b are combined to make the likelihood variables. We impose 
requirements on the likelihood ratio 7Z s /h = £si g / (£si g +£bk g ) to maximize the figure-of-merit 
(FoM) defined as A^ c / ^jN^p + N^ g , where N^p (iVy^) represents the expected number of 

signal (background) events in the signal region. We estimate A^| c using Monte Carlo (MC) 
events, while iV bkg is determined from events outside the signal region. The requirement 
for lZ s /b depends both on the decay mode and on the flavor-tagging quality, r, which is 
described in Sec. II H. The threshold values range from 0.1 (used for r > 0.875) to 0.4 (used 
for r < 0.25) for the decay B° — > <pKg (Kg — > tt+tt"), and from 0.25 to 0.65 for the decay 
B° — > K + K~K S . For the B° — »■ K + K~K S candidates, we also require | cos# T | < 0.9 prior 
to the 7£ s /b requirement. We impose a more stringent requirement, 1Z s /b > 0.75, for all r 
values in the decay B° — > (pK° s (Kg — > 7r°7r°). The 7Z s /h requirement reduces the continuum 
background by 73% for B° -> <pK s (K° s -> tt+tt"), 92% for B° -> K + K~K S and 94% for 
B° -> <PK° S (K° s -> 7T 7r°), retaining 91% of the signal for B° -> (f)K s (K° s -> tt+tt-), 72% 
for B° -> K + K-K° s and 71% for B° -> <f)K° s (K° s -> tt tt ). 

We use events outside the signal region as well as a large MC sample to study the 
background components. The dominant background is from continuum. The contributions 
from BB events are small. We estimate the contamination of B° — > K + K~Kg and B° — > 
fo(980)K° s (/ (980) -> K + K~) decays in the B° -> 0^ sample from the 5 yields in K+fs:- 
mass sideband data. The contamination of B° — > K + K~K$ events in the 5° — > (pK^ 
sample is 7.1 ± 1.6% (6.2 ± 2.0%) for DS-I (DS-II). Backgrounds from the decay B° -> 
70(980)7^5 (/o(980) — > K + K~), which has a CP eigenvalue opposite to 0-^, are found to 
be 0.4i^% (0.0±g;g%) for DS-I (DS-II). The influence of these backgrounds is treated as a 
source of systematic uncertainty. 

Figures 1(a) and (c) show the Mb c distributions for the reconstructed B° — > (pK® and 
K + K~Kg candidates within the AE signal regions after flavor tagging and vertex recon- 
struction. The AE distributions for the B° — > <pKg and K + K~K§ candidates within the 
Mbc signal region are shown in Fig. 2(a) and (b), respectively. The signal yield is determined 
from an unbinned two-dimensional maximum-likelihood fit to the AE-M\> C distribution in 
the fit region defined as M bc > 5.2 GeV/c 2 for all modes, and -0.12 GeV < AE < 0.25 GeV 
for the B° -> <pK° s (K° s -> tt+tt") or K + K-K° s decay and -0.15 GeV < AE < 0.25 GeV 
for the B° -> <pK° s (K° s -> tt tt°) decay. The <pK° s (K° s -> tt+tt") signal distribution is 
modeled with a Gaussian function (a sum of two Gaussian functions) for M bc (AE). The 
(pKg (Kg — > tt°tt ) signal distribution is modeled with a smoothed histogram obtained from 
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MC events. For the continuum background, we use the ARGUS parameterization [17] for 
M bc and a linear function for AE. The fits yield 139 ± 14 B° -> $K Q S events and 398 ± 28 
B° — ^ K + K~Kg events in the signal region, where the errors are statistical only. 

C. B° <j>K° L 

Candidate — » K + K~ decays are selected with the criteria described above. We select 
K\ candidates based on KLM and ECL information. There are two classes of K\ candidates, 
which we refer to as KLM and ECL candidates. The requirements for the KLM candidates 
are the same as those used in the B° — > J/ ipK^ selection for the sin 20! measurement [8] . 
ECL candidates are selected from ECL clusters using a K\ likelihood ratio [8], which is 
calculated from the following information: the distance between the ECL cluster and the 
closest extrapolated charged track position; the ECL cluster energy; Eg/E25, the ratio of 
energies summed in 3 x 3 and 5x5 arrays of CsI(Tl) crystals surrounding the crystal at 
the center of the shower; the ECL shower width and the invariant mass of the shower. The 
likelihood ratio is required to be greater than 0.8. For both KLM and ECL candidates, we 
also require that the cosine of the angle between the K\ direction and the direction of the 
missing momentum of the event in the laboratory frame be greater than 0.6. 

Since the energy of the K\ is not measured, M^ c and AE cannot be calculated in the 
same way that is used for the other final states. Using the four-momentum of a reconstructed 
candidate and the K\ flight direction, we calculate the momentum of the K\ candidate 
requiring AE = 0. We then calculate p c ™, the momentum of the B candidate in the cms, 
and define the B meson signal region as 0.2 GeV/c < pg" s < 0.5 GeV/c. We impose 
the requirement 1Z s /b > 0.98 to reduce the continuum background by 99.4%. The signal 
efficiency of the TZ s /b requirement is 28.1%. Here TZ s /b is based on the discriminating variables 
used for the B° — > <f)Kg decay and the number of tracks originating from the IP with a 
momentum above 0.1 GeV/c. The TZ s /b requirement is chosen to optimize the FoM, which is 
calculated taking the background from both continuum and generic B decays into account. 
The K\ detection efficiency difference between data and MC is studied using the decay 
B° — > J/ipK®, and corrections are applied to the B° — * (fiK® MC events to calculate the 
FoM. If there is more than one candidate -B — > (pK^ decay in the signal region, we take 
the one with the highest 1Z s /b value. ECL candidates are not used if there is a candidate 
B° — > (f)Kl decay with a KLM candidate. We find that about 90% of signal events are 
reconstructed with KLM candidates. 

We study the background components using a large MC sample as well as data taken 
with cms energy 60 MeV below the nominal T(4S) mass (off-resonance data). The dominant 
background is from continuum. A MC study with the efficiency correction obtained from 
B° — > J/i/jK® data yields 9±5 background events from B decays, which include B° — > <fiK*°, 
4>Kg and B + — ■> (fiK* + decays. The influence of these backgrounds, including their CP 
asymmetries, is treated as a source of systematic uncertainty. 

The p c ™ distribution after flavor tagging and vertex reconstruction is shown in Fig. 1(b). 
The signal yield is determined from an extended unbinned maximum-likelihood fit in the 
range p c ™ < 1 GeV/c. The B° — > <f)K® signal shape is obtained from MC events. Back- 
ground from BB pairs is also modeled with MC. We fix the ratio between the signal yield 
and the BB background based on known branching fractions and reconstruction efficiencies; 
the uncertainty in the ratio is treated as a source of systematic error. The continuum back- 
ground distribution is represented by a smoothed histogram obtained from MC events; we 
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confirm that the function describes the off-resonance data well. The fit yields 36 ± 15 ± 10 
B° — > §K\ events, where the first error is statistical and the second error is systematic. 
The sources of the systematic error include uncertainties in the efficiency corrections, in BB 
background branching fractions and in the background parameterizations. The result is in 
good agreement with the expected B — ► §K\ signal yield (36 ± 9 events) obtained from 
MC after applying the efficiency correction from the -B — > J/ipK L data. 

D. B° -> f (980)K° 

Candidate K s — > 7r + 7r~ decays are selected with criteria that are slightly different from 
those used for the B° — > 4>K S decay so as to obtain the best performance for the B° — > 
fo(980)Kg decay. Pairs of oppositely charged pions that have invariant masses between 
0.890 and 1.088 GeV/c 2 are used to reconstruct /o(980) — > 7r + 7r~ decays. Tracks that are 
identified as kaons (Hk/k > 0.7) or electrons are not used. We require that both K s n + and 
K s n~ combinations have invariant masses more than 0.1 GeV/c 2 above the nominal charged 
D meson mass; this removes background from — > and K #± — > K^n* 1 decays. 

The B meson signal region is defined as \AE\ < 0.06 GeV and 5.27 GeV/c 2 < M bc < 
5.29 GeV/c 2 . The AE resolution is 0.019 GeV. The dominant background is from continuum. 
For the continuum suppression, we require lZ s /b > 0.6 for events with the best-quality flavor 
tagging (r > 0.875), and 7£ s / b > 0.8 for other events. The continuum background is reduced 
by 95%, retaining 60% of signal events. Here the signal likelihood ratio TZ s /h is obtained 
from cos 6b and J 7 , which consists of the modified Fox- Wolfram moments and cos8t- 

Figure 1(d) shows the M bc distribution for the reconstructed B° — > f (980)K s candidates 
within the AE signal region after flavor tagging and vertex reconstruction. The AE distri- 
bution for the B° — > fo(980)K s candidates within the M bc signal region is shown in Fig. 2(c). 
For the signal yield extraction, we first perform an unbinned two-dimensional maximum- 
likelihood fit to the AE-M hc distribution in the fit region defined as M bc > 5.2 GeV/c 2 and 
—0.3 GeV < AE < 0.4 GeV. The signal is modeled with a Gaussian function (a sum of 
two Gaussian functions) for M bc (AE). For the continuum background, we use the AR- 
GUS parameterization for M bc and a linear function for AE. The fit yields the number of 
B° — > 7i + 7i~K s events that have n + 7i~ invariant masses within the /o(980) resonance re- 
gion, which may include contributions from B° — > p°-ft1j as well as non-resonant three-body 
B° — > it + it~K s decays. To separate these peaking backgrounds from the B° — > f (980)K s 
decay, we perform another fit to the 7t + tt~ invariant mass distribution for the events inside 
the AE-Mb c signal region. We use Breit-Wigner functions for the B° — > /q^SO)/^ signal 
as well as for B° — > pK s and a possible resonance above the / (980) mass region, which is 
referred to as /x(1300). We use the masses and widths of / (980) and /x(1300) obtained 
from data [12]. Three-body B° — > it + it~K s decays are modeled with a fourth-order poly- 
nomial function. The continuum background is modeled with a sum of a threshold function 
and a Breit-Wigner function for p resonance. The 7r + 7r _ invariant mass distribution with 
the fit result is shown in Fig. 3. The fit yields 94 ± 14 B° — > f (980)K s events. The peaking 
background contribution in the AE-M hc signal region is estimated to be 9 ± 3 events. 
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E. B° -» r)'K% 



Candidate Kg — > 7r + 7r _ decays are selected with the same criteria as those used for 
the B° — > <t>Kg decay. Charged pions from the 77, p° or 77' decay are selected from tracks 
originating from the IP. We reject kaon candidates by requiring TZk/w < 0.9. Candidate 
photons from 7r° — > 77 decays are required to have E 1 > 0.05 GeV. The reconstructed n° 
candidate is required to satisfy 0.118 GeV/c 2 < M 77 < 0.15 GeV/c 2 and p c ™ > 0.1 GeV/c, 
where M 77 and p™ s are the invariant mass and the momentum in the cms, respectively. 
Candidate photons from 77 — > 77 (77' — > p°7) decays are required to have i? 7 > 0.05 (0.1) GeV. 
The invariant mass of the photon pair is required to be between 0.5 and 0.57 GeV/c 2 for the 
rj — > 77 decay. The 7r + 7r~7r° invariant mass is required to be between 0.535 and 0.558 GeV/c 2 
for the 77 — > 7r + 7r _ 7T° decay. A kinematic fit with an 77 mass constraint is performed using 
the fitted vertex of the it + it~ tracks from the rj' as the decay point. For rj' — > p°7 decays, 
candidate p° mesons are reconstructed from pairs of vertex-constrained 7r + 7r~ tracks with 
an invariant mass between 0.55 and 0.92 GeV/c 2 . The 77' — > r]ii + 7i^ candidates are required 
to have a reconstructed mass between 0.94 and 0.97 GeV/c 2 (0.95 and 0.966 GeV/c 2 ) for 
the 77 — > 77 (77 — > 7T + 7r _ 7r ) decay. Candidate 77' — > p°7 decays are required to have a 
reconstructed mass from 0.935 to 0.975 GeV/c 2 . 

The 5 meson signal region is defined as \AE\ < 0.06 GeV for B° — > (77' — > p°7), —0.1 
GeV < A£ < 0.08 GeV for B° -> 77'Kg (77' -> rpx+Tr-, 77 -> 77), -0.08 GeV < A£ < 0.06 
GeV for B° -> 77^5 (V -> ?77r+7r-, 77 -> tt+tt-tt ), and 5.27 GeV/c 2 < M bc < 5.29 GeV/c 2 
for all decays. The AE resolution is 0.017 GeV for 77' -> p 7 , 0.027 GeV for 77' -> rpr+Tr- (77 -> 
77) and 0.018 GeV for 77' — > i]7i + ii^ (77 — > 7r + 7r _ 7r°). The continuum suppression is based 
on the likelihood ratio TZ s /b obtained from the same discriminating variables used for the 
B° — > (pKg decay, except that we only use cos9 H for the decay 77' — > p7 (p — > 7r + 7r _ ), 
where ^ is defined as the angle between the 77' meson momentum and the daughter 7r + 
momentum in the p meson rest frame. The minimum 1Z s /i, requirement depends both on 
the decay mode and on the flavor-tagging quality, and ranges from (i.e., no requirement) 
to 0.4. For the 77' — > p°7 mode, we also require | cos#t| < 0.9 prior to the TZ s /h requirement. 
With these requirements, the continuum background is reduced by 87% for 77' — > p°7, 58% 
for 77' — > 777r + 7r~ (77 — > 77) and 31% for 77' — > 777r + 7r~ (77 — > 7r + 7r~7r°), retaining 78% of the 
signal for 77' — * p°7, 94% for 77' — > 777r + 7r~ (77 — > 77) and 97% for 77' — > 777r + 7r~ (77 — > 7r + 7r _ 7r°) 

We use events outside the signal region as well as a large MC sample to study the back- 
ground components. The dominant background is from continuum. In addition, according 
to MC simulation, there is a small (~ 3%) combinatorial background from BB events in 
B° — > rj'Kg (77' — > p°7). The contributions from BB events are smaller for other modes. 
The influence of these backgrounds is treated as a source of systematic uncertainty. 

Figure 1(e) shows the M^ c distribution for the reconstructed B° — > 77'^^ candidates 
within the AE signal region after flavor tagging and vertex reconstruction, where all sub- 
decay modes are combined. The AE distribution for the B° — > rj'Kg candidates within the 
Mbc signal region is shown in Fig. 2(d). The signal yields are determined from unbinned 
two-dimensional maximum-likelihood fits to the AE-Mb c distributions in the fit region de- 
fined as M bc > 5.2 GeV/c 2 and -0.25 GeV < AE < 0.25 GeV. We perform the fit for each 
final state separately. The r/'Kg signal distribution is modeled with a sum of two (three) 
Gaussian functions for M bc (AE). For the continuum background, we use the ARGUS pa- 
rameterization for Mb c and a linear function for AE. For the 77' — > p7 mode, we include 
the BB background shape obtained from MC in the fits. The fits yield a total of 512 ± 27 
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B° — > rj'Kg events in the signal region, where the error is statistical only. 

F. B° -> ujK% 

Candidate K$ — > 7r + 7r~ decays are selected with criteria that are identical to those used 
for the B° — > 4>Kg decay. Pions for the uo — > 7r + 7r~7r° decay are selected with the same 
criteria used for the 77 — > 7r + 7r~7r° decay, except that we require p™ s > 0.35 GeV/c. The 
7T + 7r~7r° invariant mass is required to be within 0.03 GeV/c 2 of the nominal uo mass. The B 
meson signal region is defined as —0.10 GeV < AE < 0.08 GeV GeV and 5.27 GeV/c 2 < 
Mb c < 5.29 GeV/c 2 . The AE resolution is 0.028 GeV. The dominant background is from 
continuum. The continuum suppression is based on the likelihood ratio 7Z s /b obtained from 
the same discriminating variables used for the B° — > (f)Kg decay; the helicity angle $h is 
defined as the angle between the B° meson momentum and the cross product of the 7r + 
and 7T~ momenta in the uo meson rest frame. We also require | cos 8t\ < 0.9 prior to the 
7l s /b requirement. The minimum 7£ s /b requirement depends on the flavor-tagging quality, 
and ranges from 0.3 (used for r > 0.875) to 0.9 (used for r < 0.25). The lZ s /b and | cos#t| 
requirements reject 98% of the continuum background while retaining 56% of the signal. 
The contribution from BB events is negligibly small. 

Figure 1(f) shows the M bc distribution for the reconstructed B° — > 00 Kg candidates within 
the AE signal region after flavor tagging and vertex reconstruction. The AE distribution 
for the B° — > uoKg candidates within the M^ c signal region is shown in Fig. 2(e). The 
signal yield is determined from an unbinned two-dimensional maximum-likelihood fit to the 
AE-M hc distribution in the fit region defined as M bc > 5.2 GeV/c 2 and -0.12 GeV < AE < 
0.25 GeV. The signal distribution is modeled with a sum of two (three) Gaussian functions 
for Mbc {AE). For the continuum background, we use the ARGUS parameterization for 
Mb c and a linear function for AE. The fit yields 31 ± 7 B° — > uoKg events in the signal 
region with a statistical significance (£) of 7.3, where £ is defined as £ = ^—2 \n(Co/ CN aig ), 
and Cq and Cjy sis denote the maximum likelihoods of the fits without and with the signal 
component, respectively. 

G. B° -> K°7T° 

Candidate Kg — > 7r + 7r~ decays are selected with the same criteria as those used for the 
B° — > <pK s decay, except that we impose a more stringent invariant mass requirement; only 
pairs of oppositely charged pions that have an invariant mass within 0.015 GeV/c 2 of the 
nominal K° s mass are used. The 7r° selection criteria are the same as those used for the 
B° -> r]'K° s decay. 

The B meson signal region is defined as —0.15 GeV < AE < 0.1 GeV and 5.27 GeV/c 2 < 
M bc < 5.29 GeV/c 2 . The AE distribution for T^tt has a tail toward lower AE. The AE 
resolution is 0.047 GeV for the main component. The width of the tail is about 0.1 GeV. The 
dominant background is from continuum. In addition, according to MC simulation, there 
is a small (~ 2%) contamination from other charmless rare B decays. We use extended 
modified Fox- Wolfram moments, which were applied for the selection of the B° — > 7r 7r 
decay [18], to form a Fisher discriminant T . We then combine likelihoods for JF and cos 8b 
to obtain the event likelihood ratio lZ s /b for continuum suppression. 
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As described below, we include events that do not have B decay vertex information in 
our fit to obtain better sensitivity for the CP- violation parameter A. For events with vertex 
information, the high-7?. s /b region is defined as 7£ s /b > 0.78 and the low-7£ s /b region as 
0.4 < 1Z s / h < 0.78. For events without vertex information, the high-7?. s /b region is defined 
as TZ s/h > 0.74 (0.76) for DS-I (DS-II), and the low-TZ s/h region as 0.4 < TZ s/h < 0.74 (0.76) 
for DS-I (DS-II). By the high-7?. s /b requirement, 95% of continuum backgrounds are rejected 
and 64% of signal events remain. As for the low-7£ s /b region, 84% of continuum backgrounds 
are rejected and 26% of signal events remain. 

Figure 1(g) shows the M bc distribution for the high-7?. s /b B° — > Kgir candidates within 
the AE signal region after flavor tagging and before vertex reconstruction. Also shown 
in Fig. 1(h) is the M^ c distribution for the low-7£ s /b B° — > Kgir candidates. The AE 
distributions for the high-7?. s /b and low-7£ s /b candidates are shown in Fig. 2(f) and (g), 
respectively. The signal yield is determined from an unbinned two-dimensional maximum- 
likelihood fit to the AE-M^ C distribution in the fit region defined as 5.2 GeV/c 2 < M^ c < 
5.29 GeV/c 2 and -0.5 GeV < AE < 0.5 GeV. The B° -> K° s ti signal distribution is 
modeled with a Gaussian function for M bc and with a Crystal Ball function [19] for AE. 
For the continuum background, we use the ARGUS parameterization for M bc and a second- 
order Chebyshev function for AE. The B decay background distribution is represented by 
a smoothed histogram obtained from MC simulation. The fits yield 168 ± 16 and 79 ± 19 
B° — > Kgii events in the high-7?. s /b and low-7?. s /b signal regions, respectively, where the 
errors are statistical only. The same procedure after the vertex reconstruction yields a total 
of 71 ± 13 K s n° events. 

H. Flavor Tagging 

The 6-flavor of the accompanying B meson is identified from inclusive properties of par- 
ticles that are not associated with the reconstructed B° — > f C p decay. We use the same 
procedure that is used for the sin 20! measurement [9]. The algorithm for flavor tagging is 
described in detail elsewhere [20]. We use two parameters, q and r, to represent the tagging 
information. The first, q, is defined in Eq. (1). The parameter r is an event-by-event, MC- 
determined flavor-tagging dilution factor that ranges from r = for no flavor discrimination 
to r = 1 for unambiguous flavor assignment. It is used only to sort data into six r intervals 
listed in Table I. The wrong tag fractions for the six r intervals, wi (I = 1,6), and differ- 
ences between B° and B° decays, Awi, are determined from the data; we use the same values 
that were used for the sin20i measurement [9] for DS-I. Wrong tag fractions for DS-II are 
separately obtained with the same procedure and are listed in Table I. The total effective 
tagging efficiency for DS-II is determined to be e e g = X)f=i e^(l — 2w/) 2 = 0.30 ± 0.01, where 
ei is the event fraction for each r interval determined from the J/ipK s data and is listed in 
Table I. The error includes both statistical and systematic uncertainties. We find that the 
wrong tag fractions for DS-II are slightly smaller than those for DS-I. As a result, the e e gp 
value for DS-II is slightly larger than that for DS-I (e eff = 0.287 ± 0.005). 

I. Vertex Reconstruction 

The vertex position for the f C p decay is reconstructed using charged tracks that have 
enough SVD hits; at least one layer with hits on both sides and at least one additional z 
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hit in other layers for SVD-I, and at least two layers with hits on both sides for SVD-II. A 
constraint on the IP is also used with the selected tracks; the IP profile is convolved with 
the finite B flight length in the plane perpendicular to the z axis. The pions from K s decays 
are not used except in the analysis of B° — > K s n° decays. The typical vertex reconstruction 
efficiency and z resolution for B° — > 4>K S decays are 95% and 78 /iin, respectively. Similar 
values are obtained for other fcp decays except for B° — > K s tt° decays. 

The vertex for B° — > K s n° decays is reconstructed using the K s trajectory and the IP 
constraint, where both pions from the K s decay are required to have enough SVD hits in 
the same way as that for other f C p decays. The reconstruction efficiency depends both on 
the K s momentum and on the SVD geometry; the efficiency with SVD-II (32%) is higher 
than that with SVD-I (23%) because of the larger outer radius and the additional layer. 
The typical z resolution of the vertex reconstructed with the K s is 93 /xm for SVD-I and 
110 jum for SVD-II. 

The /tag vertex determination with SVD-I remains unchanged from the previous pub- 
lication [4], and is described in detail elsewhere [21]; to minimize the effect of long-lived 
particles, secondary vertices from charmed hadrons and a small fraction of poorly recon- 
structed tracks, we adopt an iterative procedure in which the track that gives the largest 
contribution to the vertex x 2 is removed at each step until a good x 2 is obtained. The 
reconstruction efficiency was measured to be 93%. The typical z resolution is 140/im [8]. 

For SVD-II, we find that the same vertex reconstruction algorithm results in a larger 
outlier fraction when only one track remains after the iteration procedure. Therefore, in 
this case, we repeat the iteration procedure with a more stringent requirement on the SVD- 
II hit pattern; at least two of the three outer layers have hits on both sides. The resulting 
outlier fraction, which is described in Sec. Ill, is comparable to that for SVD-I, while the 
inefficiency caused by this change is small (2.5%). 

J. Summary of Signal Yields 

The signal yields for B° — > f CP decays, 7V sig , after flavor tagging and vertex reconstruction 
(before the vertex reconstruction for the decay B° — > K s tt°) are summarized in Table II. 
The signal purities are also listed in the table. 

III. RESULTS OF CP ASYMMETRY MEASUREMENTS 

We determine S and A for each mode by performing an unbinned maximum-likelihood 
fit to the observed At distribution. The probability density function (PDF) expected for 
the signal distribution, V S i g (At; S, A, q, wi, Awi), is given by Eq. (1) incorporating the effect 
of incorrect flavor assignment. The distribution is convolved with the proper-time interval 
resolution function R sig (At), which takes into account the finite vertex resolution. 

For the decays B° -> <j>K%, K + K~K S , <f)K L , f {980)K s , rj'K s and uK s , we use flavor- 
specific B decays governed by semileptonic or hadronic b — > c transitions to determine the 
resolution function. We perform a simultaneous multiparameter fit to these high-statistics 
control samples to obtain the resolution function parameters, wrong-tag fractions (Sec- 
tion II H), Arrid, t b + and t b o. We use the same resolution function used for the sin 20! 
measurement for DS-I [9]. For DS-II, the following modifications are introduced: a sum 
of two Gaussian functions is used to model the resolution of the fcp vertex while a single 
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Gaussian function is used for DS-I; a sum of two Gaussian functions is used to model the 
resolution of the tag-side vertex obtained with one track and the IP constraint, while a 
single Gaussian function is used for DS-I. These modifications are needed to account for dif- 
ferences between SVD-I and SVD-II, as well as different background conditions in DS-I and 
DS-II. We test the new resolution parameterization using MC events on which we overlay 
beam-related background taken from data. A fit to the MC sample yields correct values for 
all parameters. With the multiparameter fit to data, we find that the standard deviation of 
the main Gaussian component of the resolution function is reduced from 78 jum to 55 /xm, 
which is consistent with our expectation from the improved impact parameter resolution of 
SVD-II [15]. The same fit also yields r B « = 1.519 ± 0.010 ps, t b+ = 1.652 ± 0.012 ps and 
Arrid = 0.516 ± 0.007 ps -1 , where errors are statistical only. The results are consistent with 
those obtained with DS-I [9] and also with the world average values [22]. Thus we conclude 
that the resolution of SVD-II is well understood. 

For the decay B° -> K%n°, we use the resolution function described above with additional 
parameters that rescale vertex errors. The rescaling parameters depend on the detector 
configuration (SVD-I or SVD-II), SVD hit patterns of charged pions from the Kg decay, and 
K s decay vertex position in the plane perpendicular to the beam axis. These parameters 
are determined from a fit to the At distribution of B° — > J/^K S data. Here the Kg and 
the IP constraint are used for the vertex reconstruction, the B° lifetime is fixed at the world 
average value, and 6-flavor tagging information is not used so that the expected PDF is an 
exponential function convolved with the resolution function. 

We check the resulting resolution function by also reconstructing the vertex with leptons 
from J/ip decays and the IP constraint. We find that the distribution of the distance 
between the vertex positions obtained with the two methods is well represented by the 
obtained resolution function convolved with the well-known resolution for the J/ijj vertex. 
Finally, we also perform a fit to the B° — > J/ipK s sample with 6-flavor information and 
obtain Sj/^ K o = +0.68 ± O.lO(stat) and Aj/^ K o = +0.02 ± 0.04(stat), which are in good 
agreement with the world average values. Thus, we conclude that the vertex resolution for 
the B° — > K S 7T° decay is well understood. 

We determine the following likelihood value for each event: 

Pi = (l-/oi) /[/ sig nig(AOi?si g (A^-AO 

+ (1 - /si g )n kg (At , ) J Rbk g (At, - Af)]d(AO 

+ / iP i(A^), (2) 

where P Q \(At) is a broad Gaussian function that represents an outlier component with a 
small fraction f Q \ [9]. The width of the outlier component for DS-II is determined to be 
(44 ±5) ps; the fractions of the outlier components are (3.1 ± 1.2) x 10~ 4 for events with the 
/tag vertex reconstructed with more than one track, and (1.2 ±0.1) x 10 -2 for the case only 
one track is used. These values are comparable to those for DS-I [9]. The signal probability 
/si g depends on the r region and is calculated on an event-by-event basis as a function of p c ™ 
for the B° — > (frK® decay and as a function of AE and Mb c for the other modes. A PDF for 
background events, Pbkg(At), is modeled as a sum of exponential and prompt components, 
and is convolved with a sum of two Gaussians i?bk g - Parameters in P bkg (At) and i? bkg for 
continuum background are determined by the fit to the At distribution for events outside 
the AE-M hc signal region except for the B° — > §K\ decay. For the B° — > 4>K^ decay, we 
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use p c ™ sideband events to obtain the parameters. Parameters in Pbkg(At) and i?bkg for 
BB background events in B° — > rj'Kg, B° — > Kgir and B° — > §K\ decays are determined 
from MC simulation. 

We fix r B o and Am d at their world average values [22]. We assume no CP asymmetry in 
the background At distributions and possible CP asymmetries in the B decay backgrounds 
are treated as sources of systematic error. In order to reduce the statistical error on A, we 
include events without vertex information in the analysis of B° — > K^n . The likelihood 
value in this case is obtained by integrating Eq. (2) over Atj. 

The only free parameters in the final fit are S and A, which are determined by maximizing 
the likelihood function L = ]li Pi(Atf, S, A) where the product is over all events. Table III 
summarizes the fit results of S and A. We define the raw asymmetry in each At bin by 
(N q=+ i — N q= _i) I (N q=+1 + N q= _i), where iV 9=+1 (_ 1 ) is the number of observed candidates 
with q = +1(— 1) [23]. Figures 4(a-f) show the raw asymmetries in two regions of the flavor- 
tagging parameter r. While the numbers of signal events in the two regions are similar, 
the effective tagging efficiency is much larger and the background dilution is smaller in the 
region 0.5 < r < 1.0. Note that these projections onto the At axis do not take into account 
event-by-event information (such as the signal fraction, the wrong tag fraction and the vertex 
resolution), which is used in the unbinned maximum-likelihood fit. 

Tables IV and V list the systematic errors on S and A, respectively. The total systematic 
errors are obtained by adding each contribution in quadrature, and are much smaller than 
the statistical errors for all modes. 

To determine the systematic error that arises from uncertainties in the vertex recon- 
struction, the track and vertex selection criteria are varied to search for possible systematic 
biases. Small biases in the Az measurement are observed in e + e~ — > and other con- 

trol samples. Systematic errors are estimated by applying special correction functions to 
account for the observed biases, repeating the fit, and comparing the obtained values with 
the nominal results. The systematic error due to the IP constraint in the vertex reconstruc- 
tion is estimated by varying (±10/xm) the smearing used to account for the B flight length. 
Systematic errors due to imperfect SVD alignment are determined from MC samples that 
have artificial misalignment effects to reproduce impact-parameter resolutions observed in 
data. 

Systematic errors due to uncertainties in the wrong tag fractions are studied by varying 
the wrong tag fraction individually for each r region. Systematic errors due to uncertainties 
in the resolution function are also estimated by varying each resolution parameter obtained 
from data (MC) by ±la (±2<r), repeating the fit and adding each variation in quadrature. 
Each physics parameter such as r B o and Am^ is also varied by its error. A possible fit bias 
is examined by fitting a large number of MC events. 

Systematic errors from uncertainties in the background fractions and in the background 
At shape are estimated by varying each background parameter obtained from data (MC) 
by ±1(7 (±2(j). 

Additional sources of systematic errors are considered for B decay backgrounds that are 
neglected in the PDF. We consider uncertainties both in their fractions and CP asymmetries; 
for modes that have non-vanishing CP asymmetries, we conservatively vary the CP-violation 
parameters within the physical region and take the largest variation as the systematic error. 
The effect of backgrounds from K + K~K° S and f {980)K° s (/ (980) -> K+K~) in the B° -> 
4>Kg sample is considered. Uncertainties from B — > <pK* and other rare B decay backgrounds 
in the B° — > (pK L sample are also taken into account. Effects of possible CP asymmetries 
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in B decay backgrounds for 7^71"° and /o(980)-Kg are evaluated. The peaking background 
fraction in the B° — > /o(980)i£s sample depends on the functions used in the fit to the 
7r + 7r~ invariant mass distribution. The systematic errors due to the uncertainties of the 
masses and widths of the resonances used in the fit are estimated. The width of / (980) 
is varied between 0.04 GeV/c 2 and 0.1 GeV/c 2 . The mass and width of the /x(1300) are 
varied by ±0.03 GeV/c 2 and ±0.1 GeV/c 2 , respectively. The effect of possible interference 
between resonant and non-resonant amplitudes, which is neglected in the nominal analysis, 
is also evaluated. We perform a fit to the tc + it~ distribution of a MC sample generated 
with interfering amplitudes and phases for B — > Knn decays measured from data [12]. 
The observed difference in the signal yield from the true value is taken into account in the 
systematic error determination. We also repeat the fit to the At distribution ignoring the 
contribution of the peaking background. The differences in S and A from our nominal 
results are included in the systematic error. 

Finally, we investigate the effects of interference between CKM-favored and CKM- 
suppressed B — > D transitions in the /tag final state [24]. A small correction to the PDF for 
the signal distribution arises from the interference. We estimate the size of the correction 
using the B° — > D*~l + v sample. We then generate MC pseudoexperiments and make an 
ensemble test to obtain systematic biases in S and A. In general, we find effects on S are 
negligibly small, while there are sizable possible shifts in A. 

Various crosschecks of the measurements are performed. We reconstruct charged B meson 
decays that are the counterparts of the B° —> f C p decays and apply the same fit procedure. 
All results for the S term are consistent with no CP asymmetry, as expected. Lifetime 
measurements are also performed for the fcp modes and the corresponding charged B 
decay modes. The fits yield r B a and r B + values consistent with the world average values. 
MC pseudoexperiments are generated for each decay mode to perform ensemble tests. We 
find that the statistical errors obtained in our measurements are all consistent with the 
expectations from the ensemble tests. 

For the B° -> <pK decay, a fit to DS-I alone yields S = -0.68±0.46(stat) and A = 0.00 ± 
0.28(stat), while a fit to DS-II alone yields S = ±0.80±0.45(stat) and A = ±0.15±0.33(stat). 
Note that the results for DS-I differ from our previously published results S = —0.96 ± 
0.501^? and A = -0.15 ± 0.29 ± 0.07 [4], as the decays B° -> <f)K L and <j>K° s (K° s -> 7r°7r°) 
are included in this analysis. Fit results to B° <pK s (K s — > 7r + 7r~) decays in DS-I in 
this analysis are S = —0.97 ± 0.50(stat) and A — 0.15 ± 0.29(stat), which are consistent 
with the previous results. From MC pseudoexperiments, the probability that the difference 
between the S values in DS-I and DS-II is larger than the observed difference (1.46) is 
estimated to be 4.5%. Fits to only <pK s decays and only <pK L decays yield S = 0.02 ± 0.33 
and A = 0.07 ± 0.22, and -S = 2.3 ± 2.0 and A = 0.6 ± 1.2, respectively. A sin20i 
measurement with DS-II is performed using B° — > J/ipK s (K s — > it + it~ or 7r°7r°) and 
B° — > J/ipK L decays as a crosscheck. Applying the same procedure to both DS-I and DS-II, 
we obtain Sj/^ K o = ±0.696 ± 0.061(stat) and Aj/^ K o = ±0.011 ± 0.043(stat) for DS-I, and 
Sj/^ K o = +0.629 ±0.069(stat) and A J/i)K o = ±0.035 ±0.044(stat) for DS-II. The results are 
in good agreement with each other, and are also consistent with SM expectations. As all the 
other checks mentioned above also yield results consistent with expectations, we conclude 
that the difference in S^k between the two datasets is due to a statistical fluctuation. 

For the B° — > i]'K s , the statistical significance of the CP asymmetry is calculated with 
the Feldman- Cousins frequentist approach [25]. The case with no CP asymmetry (5 = 
and A = 0) is ruled out at 99.92% confidence level, equivalent to 3.4 standard deviations 
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for Gaussian errors. 

As discussed in Section I, to a good approximation, the SM predicts S = — £/sin20i 
for the B° -> <f>K°, f {980)K° s , rfK° s , uK° s and K° s n° decays. For the B° -> K + K-R S 
decay, the SM prediction is given by S = — (2/ + — 1) sin 20!. The effective sin 20! for this 
mode is obtained to be 0.74 ± 0.27 ± 0.06±g;f|. The third error is an additional systematic 
error arising from the uncertainty of the = +1 fraction. Figure 5 summarizes the sin 20! 
determination based on our S measurements for these decays. For each mode, the first 
error shown in the figure is statistical and the second error is systematic. We also include 
the result of the time-dependent CP- violating asymmetry measurement in B° — > KgKgKg 
decays by Belle [26]. We obtain sin 20i = +0.40±0.13 as a weighted average, where the error 
includes both statistical and systematic errors. The result differs from the SM expectation 
by 2.4 standard deviations. 

IV. SUMMARY 

We have performed improved measurements of CP-violation parameters for B° — ► (f)K° 
(including both 4>Kg and 0A^), K + K~K$ and 7]' Kg decays, and new measurements for 
B° — > f (980)Kg, ujK s and A^7r° decays. These charmless decays are dominated by b — > s 
flavor- changing neutral currents and are sensitive to possible new CP- violating phases. The 
results for each individual decay mode are consistent with the SM expectation within two 
standard deviations except for the B° — > f (980)Kg decay. The combined result for the 
B° -> <pK°, K + K~K° S , f (980)K° s , r]'K° s , uK° s , K° s ir° and previously reported K° S K° S K° S 
decays differs from the SM expectation by 2.4 standard deviations. Measurements with 
a much larger data sample are required to conclusively establish the existence of a new 
CP-violating phase beyond the SM. 
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TABLE I: Event fractions e;, wrong-tag fractions Wi, wrong-tag fraction differences Awi, and 
average effective tagging efficiencies e 1 ^ = e/(l — 2wi) 2 for each r interval for DS-II. Errors for wi 
and Awi include both statistical and systematic uncertainties. The event fractions are obtained 
from J/ipK® data. 



1 


r interval 


e/ 


Wl 


Awi 


e eff 


1 


0.000 


- 0.250 


0.397 ±0.015 


0.464 ± 0.007 


±0.008 ± 0.007 


0.002 ± 0.001 


2 


0.250 


- 0.500 


0.146 ±0.009 


0.320 ± 0.008 


-0.022 ±0.010 


0.019 ±0.002 


3 


0.500 


- 0.625 


0.108 ±0.008 


0.225 ±0.011 


±0.029 ±0.011 


0.033 ± 0.004 


4 


0.625 


- 0.750 


0.107 ±0.008 


0.158 ±0.010 


±0.003 ±0.011 


0.050 ± 0.005 


5 


0.750 


- 0.875 


0.098 ± 0.007 


0.109 ±0.009 


-0.028 ±0.011 


0.060 ± 0.005 


6 


0.875 


- 1.000 


0.144 ±0.009 


0.015 ±0.005 


±0.007 ± 0.007 


0.135 ±0.009 



TABLE II: Estimated signal purities and signal yields 7V sig in the signal region for each fcp mode 
that is used to measure CP asymmetries. The result for the B° — > K^ir decay is obtained with 
the sample after flavor tagging but before vertex reconstruction. The results for other decays are 
obtained after flavor tagging and vertex reconstruction. 



Mode 




purity A^sig 


<j>Kl 

K+K-K° s 
/ (980)K° 
ri'K° s 
uK° s 


(high-7e s/b ) 
(low-7e s/b ) 


-1 
+1 

+1 (83%), -1 (17%) 

+ 1 

-1 

-1 

-1 

-1 


0.63 139 ± 14 
0.17 36 ±15 
0.56 398 ± 28 
0.53 94 ± 14 
0.61 512 ±27 
0.56 31 ±7 
0.54 168 ± 16 
0.15 79 ±19 


TABLE III: Results of the fits to the At distributions. The first error is statistical and the second 
error is systematic. We combine B° — > 4>K^ and B° — > 4>K\ decays to obtain S^ K o and A^o. 


Mode 


SM expectation for S 


S 


A 


K + K-K Q S 
/ (980)K° 

7]'K° S 


±sin 20i 
-(2/+-l)sin20 1 
—sin 20i 
+sin 20i 
+sin 20i 
+sin20i 


+0.08 ±0.33 ±0.09 
-0.49 ±0.18 ±0.04 
+0.47 ±0.41 ±0.08 
+0.65 ±0.18 ±0.04 

+0.76 ± 0.65lJu6 
+0.32 ±0.61 ±0.13 


+0.08 ±0.22 ±0.09 
-0.09 ±0.12 ±0.07 
-0.39 ±0.27 ±0.09 
-0.19 ±0.11 ±0.05 
+0.27 ±0.48 ±0.15 
-0.11 ±0.20 ±0.09 
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TABLE IV: Summary of the systematic errors on S. 





(j)K v 


K + K 


K s 


/o(980)Xg 


rj'K° s 


uK° s 


t^o o 


Vertex reconstruction 


0.01 




0.01 


0.02 


0.01 


0.01 


0.02 


Flavor tagging 


0.01 


< 


0.01 


0.01 


0.01 


0.04 


0.01 


Resolution function 


0.04 




0.03 


0.03 


0.03 


0.07 


0.05 


Physics parameters 


< 0.01 


< 


0.01 


0.01 


< 0.01 


0.01 


0.02 


Possible fit bias 


0.01 




0.01 


0.03 


0.01 


+0.01 
-0.10 


0.02 


Background fraction 


0.08 




0.02 


0.05 


0.02 


0.10 


0.08 


Background At shape 


0.01 


< 


0.01 


0.04 


< 0.01 


0.02 


0.08 


Tag-side interference 


< 0.01 


< 


0.01 


< 0.01 


< 0.01 


0.01 


< 0.01 


Total 


0.09 




0.04 


0.08 


0.04 


+0.13 
-0.16 


0.13 



TABLE V: Summary of the systematic errors on A. 





4>K° 


K+K 


-K° 


/o(980)K° 




uK% 


K^ 


Vertex reconstruction 


0.03 




0.04 


0.04 


0.03 


0.04 


0.04 


Flavor tagging 


< 0.01 


< 


0.01 


0.01 


0.01 


< 0.01 


0.01 


Resolution function 


0.02 




0.02 


0.02 


0.01 


0.04 


< 0.01 


Physics parameters 


< 0.01 


< 


0.01 


< 0.01 


< 0.01 


< 0.01 


< 0.01 


Possible fit bias 


0.01 




0.01 


0.03 


0.01 


+0.01 
-0.03 


0.01 


Background fraction 


0.04 




0.01 


0.06 


0.02 


0.14 


0.04 


Background At shape 


0.03 


< 


0.01 


0.01 


< 0.01 


0.03 


0.05 


Tag-side interference 


0.06 




0.06 


0.04 


0.03 


0.03 


0.05 


Total 


0.09 




0.07 


0.09 


0.05 


0.15 


0.09 
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FIG. 1: M bc distributions for (a) B° -» <j>K%, (c) B° -► K + K~ K% (d) 5° -» / (980)A% (e) 



5° -» j/Jtg, (f) 5° -► wiTg, (g) B° -► A^tt (high-ft s/b ), and (h) 5° -► A^rr (low-7£ s/b ), within 
the AA signal region and (b) p™ s distribution for B — > Solid curves show the fits to signal 

plus background distributions, and dashed curves show the background contributions. 
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FIG. 2: AE distributions for (a) B° -> <pK s , (b) B° -> K + K~K° S , (c) 5° -» f (980)K%, (d) 
5° -» r?'^, (e) B° -» cjETg, (f) 5° -» Kgvr (high-7e s/b ), and (g) 5° -> K^tt (low-7£ s/b ), within 
the Mbc signal region. Solid curves show the fit to signal plus background distributions, and dashed 
curves show the background contributions. 
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FIG. 3: Invariant mass distribution of ir + ir pairs for B° -» / (980)ifg candidates in the AE-M hc 
signal region. The solid curve shows the fit result. The dashed line shows continuum background. 
The dotted line shows the sum of the continuum and the non-resonant B° — > ir + ir^K^. The dot- 
dashed lines show the B° -> /j ^ (left), the 5° -» / (980)Kg (middle), and the 5° -> / x (1300)lTg 
(right). 
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FIG. 4: Raw asymmetry in each At bin with < r < 0.5 (top) and with 0.5 < r < 1.0 (bottom) 
for (a) B° -» cpK°, (b) B° -» K + K~ Kg, (c) 5° -» / (980)^, (d) 5° -► r/Kg, (e) B° -> cuKg, 
and (f) — > -f^^vr . The solid curves show the results of the unbinned maximum-likelihood fits. 
The dashed curves show the SM expectation with sin20i = +0.73 and .4 = 0. 
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Belle Results on Time-dependent CP Violation in b— >s 
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sin^WA) 
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_L 
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FIG. 5: Summary of sin2^>i measurements performed with hadronic B decays governed by the 
b — » sqq transition. The world average sin2^i value obtained from B° — > J/ipK° and other related 
decay modes governed by the b — > ccs transition [11] is also shown as the SM reference. 
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